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Abstract 

We present a direct and simple method for the computation of the total scat- 
tering matrix of an arbitrary finite noncompact connected quantum graph given its 
metric structure and local scattering data at each vertex. The method is inspired 
by the formalism of Reflection- Transmission algebras and quantum field theory 
on graphs though the results hold independently of this formalism. It yields a 
simple and direct algebraic derivation of the formula for the total scattering and 
has a number of advantages compared to existing recursive methods. The case 
of loops (or tadpoles) is easily incorporated in our method. This provides an 
extension of recent similar results obtained in a completely different way in the 
context of abstract graph theory. It also allows us to discuss briefly the inverse 
scattering problem in the presence of loops using an explicit example to show that 
the solution is not unique in general. On top of being conceptually very easy, the 
computational advantage of the method is illustrated on two examples of "three- 
dimensional" graphs (tetrahedron and cube) for which other methods are rather 
heavy or even impractical. 
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1 Introduction 



Excitement in the study of systems on quantum graphs has been revived recently as 
they provide models for the study of transport properties in quantum wires connected 
through junctions. It is largely motivated by the range of different physical applications 
that can be linked to such models, starting from condensed matter experiments or atomic 
wires up to chaos and neural networks, for reviews, see e.g. [U [2J. 

A powerful formalism in this respect is that of quantum fields theory on graphs 
combined with bosonization techniques. One of the central objects in this approach 
is total scattering matrix of the graph and the knowledge of its analytic structure. A 
number of results is already available in [3], SI El E] but essentially for star graphs. Results 
that apply to more general graphs can be found in [71 El El QUI Ell E21 E32]- However, all 
the different techniques presented in those papers imply the use of a recursive approach 
that becomes rather cumbersome when the graph is complicated. 

The goal of this paper is to provide an efficient and simple techniques to compute this 
matrix for an arbitrary finite noncompact connected quantum graph knowing only its 
metric structure and local scattering data at each vertex. The point of view taken here is 
that the complete graph is obtained by assembling star graphs (single vertex graphs with 
a certain number of edges) which are well-understood. We obtain an explicit formula 
for the total scattering matrix. It turns out that our results hold beyond the context 
of quantum field theory on graphs. Not only do they represent an extension of recent 
results [TT] to the case of graphs with loops but our method also provides a direct (as 
opposed to recursive [101 and simple algebraic derivation. 

The paper is organised as follows. In section [2] we pesent our formalism to compute 
directly the scattering matrix associated to a general quantum graph without loop. Once 
the notation is settled, the calculation is very simple and effective. In the next section, 
we show how to extend the techniques to graphs with loops. Then, in section m we 
illustrate the techniques in computing the scattering matrix for graphs corresponding to 
Platonic solids, the cases of tetrahedron and cube being treated in great details. Finally, 
the last section is devoted to a short conclusion on possible applications. 

2 General setting and results 

We consider a finite noncompact graph with N vertices that we label with a — 1, . . . , N 
and with internal and external edges. The graph is compact if it has no external edges. 
At each vertex a are attached possibly several edges. One can endow the graph with a 
metric structure: the external edges are associated to infinite half-lines and are connected 
to a unique vertex; the internal edges are associated to intervals of finite length and 
connect two vertices, possibly not distinct. In the case where an internal edge connects 
the same vertex, we call it a loop (also called tadpole in the literature). Two edges are 
adjacent if they are connected by an internal edge. We consider a connected graph i.e. 
a graph such that for any two vertices a, f3 there is a sequence {a,\ = a, a 2 , ■ ■ ■ , a q = (3} 
of adjacent vertices. We define an orientation on the edges, and in the case of internal 
edges, (a/3) will define an edge with orientation from vertex a to vertex (3. By convention, 
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external edges (aO) are always oriented from the vertex to infinity. On each of these 
edges, we attach modes (of fields living on the edge) 



of(p) j = l,...,N af3 ; P = 0,1,..., N ; a = l,..., N ; a ^(3, 



p being an orientation dependent parameter which has the interpretation of a momen- 
tum or a rapidity in applications to quantum fields on graphs and with the following 
conventions: 

• a = 1,2, . . . , N denotes the vertex to which the edge is attached; 

• (3 = 0, 1, 2, . . . , N denotes the vertex linked to a by the edge under consideration, 
with the convention that external edges corresponds to (3 = 0; 

• j = 1, . . . , N a p numbers the different edges between a and (3, N a p being their total 
number. We set N a/3 = if a is not connected to (3. 

In this way the ordered triplet (a,(3,j) uniquely defines all the oriented edges of the 
graph. Obviously, (a,/3,j) and ((3,a,j) define the same edge, but with a different 
orientation. Hence we have N a p = Np a . We will call internal mode (resp. external 
mode) a mode living on an internal edge (resp. external edge). 

2.1 General case without loops 

For the time being, we assume N aa = for all a = 1,. . . ,N i.e. we do not consider 
loops. We will see later on that they are easily incorporated in our formalism. The modes 
are not independent but are related by two types of fundamental relations defining the 
scattering and propagation on the graph: 

• Local scattering at vertex a: Following the RT-algebra formalism (see e.g. [T^[T5]). 



a f (P) = ££*?» < 7 (-P) Vj = 1, . . . , N aP ; V/3 = 0, 1, . . . , N (2.1) 



where s^ k (p) are the components of the local scattering matrix S a (p) which sat- 
isfies S a (p)S a (—p) = 2. 

• Propagation on edge {at(3j)\ As already mentionned, the edges {otf3j) and ([3a. j) are 



1 The particular form of this relation comes from the fact that we have in mind applications to 
quantum fields which are Fourier transforms of the modes considered here, see |10j for instance. 



this reads 



N N t 



7 =0 fc=l 
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The aim now is to obtain the scattering relations directly between the external modes 
i. e. relations of the form 



N 



,q0 



(?) 



N l0 

EE s S»«i (-p) v J = i,...,iv a0 



Va 



N. 



(2.3) 



7=1 fc=l 



where s^ t .- k {p) are the components of the total scattering matrix for the graph, Stotip)- 
This is most easily achieved by arranging the modes in vectors and using simple linear 
algebra. Denote M rX s the vector space of r x s matrices over C. In particular, we 
identify M n xn and End(C n ). We denote E\'* the r x s matrix whose only nonzero 
entry is 1 at position The set {£ , [/}i=i,..., r; is a basis of M rX s- We will drop the 

j=i,-,s 

superscripts every time this does not cause confusion i.e. each time the size of the matrix 
corresponds to the range of the indices. Similarly, we denote {e" }_,-=!,.. . jTl the canonical 
basis of C n and we will use a similar convention. Finally, we denote J-{p) the space of all 
(possibly generalized) functions of p G C, with the understanding that these functions 
can be operator- valued in quantum field applications (cf the modes). The following 
definitions illustrate our notations and conventions. For a given vertex a, we define 
different vectors: 

• We collect the external modes attached to a in 

/ af(p) 



N„ 



= e i® a f(p) e cNa0 ® Hp) 
\ <Jp) 1 3=1 



(2.4) 



• We collect the internal modes attached to a in 

/ <{P) \ 



af( P ) 



B a {p) 



a2 



(P) 



<(P) 

V <' N (p) ) 



(2.5) 
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where only the modes with N a/ 3 ^ appear. For conciseness^], we write this as 

N N a/ 3 

B a (p) =J2J2 e ?® e l® a f (?) G ^ ® ^(p) (2.6) 

0=1 j=l 

N 

where v a = £^ N a p is the number of internal edges attached to a. This makes the 
&=i 

following computations a lot more transparent but the reader should remember 
the actual content and size of each vector. 

• Similarly, we collect all the modes attached to a in 

Mp) = J2J2 e e+i ® e > ® a f (p) e £Na ® ^(p) ( 2 - 7 ) 

/3=0 j=l 

where iVa = iV Q o + is the total number of edges attached to a. This way, A a is 
the concatenation of A a and B a with A a " sitting on top" . 

With the same conventions, we introduce 

N N a0 N al 

S <*(p) = EEE ^+1,7+1 ® % ® «Sfc(p) e End(C^) ® Hv) , (2.8) 

/3,7=0 j=l k=l 

so the relations (12. ip read 

A a (p) = S a (p) A a (-p) , Va = 1, . . . , N . (2.9) 
The set of relations (I2.9p can be gathered into a single one: 

N N 

A{p) = S{p)A{-p) with A(p) = ^e a ® A a (p) and S(p) = £ E aa ® S a (p) 

a=l a=l 

(2.10) 

Remark that A{p) G C Ne+2Nl ®!F{p) where N e = £^ N a0 is the total number of external 

a=l 

edges and iVj = £^ N a/ 3 is the total number of internal edges. Then, we introduce 

l<a<P<N 

N 

B{p) = J2 e a® G £2Nl ® ^(P) > ( 2 - U ) 

a=l 



2 The explicit, longer formula is 

q a Pp+l-l N a p 
p=0 /3=/3 p + l 3=1 

where {/3i, . . . , /3 9a } are the labels /? such that N a p = and we have set (3q = and (3q a +i = N + 1 for 
convenience. 
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and 

E(p) = J2 Ea 'P ® E ^ a ® E H ® exp(-idfp) e £nd(C 2JVl ) ® ^(p) , (2.12) 

a,/3=l i=l 

so the relations (12. 2p read 

B(p) = E(p)B(-p) (2.13) 

It is easy to see that 

E(p) E(-p) = E <*>°> ® E P>P ® ^ 

that acts as the identity matrix H^at*- The final step is to decompose the matrix S(p) 
into four submatrices related to external or internal edges: 

N N a0 

S (U \p) = Y,J2 E ™® E j*® s Zk(p) eEnd(C N °)®F(p) (2.14) 

a=l j,k=l 

N N a0 Ncr, 

s {12) (p) = E ™ ® E ^ ® ^ ® s S*(p) G ^ 

N e x2Ni ® -^(p) (2-15) 

a,7=l j=l fc=l 

AT AT a0 

2NiXN e <8> -^(p) (2-16) 

a,/3=l j=l fe=l 
AT iVtt/3 V a7 

S (22) (P) = ]T ei?nrf(C 27V 0®^(p)-(2.17) 

a,/3,7=l j=l fc=l 

Therefore, the set of all the relations we have becomes 

A(p) = S {11 \p)A(-p) + S {12) (p)B(-p) (2.18) 
B(p) = S {21 \p)A(-p) + S {22 \p)B(-p) (2.19) 
B(p) = E(p)B(-p) (2.20) 

Assuming that E(j>) — S^ 22 \p) is invertible this yields the desired relations in the form 

A(p) = S tot (p)A(-p) , (2.21) 

with 

SUP) = S (11) (p) + S (12) (p) [E(p) - S^ip)]' 1 S<M(p) . (2.22) 
The internal modes can be expressed in terms of the external ones: 

B(p) = [E(-p) - S^i-p)]- 1 S (21 \-p)A(p) (2.23) 

These two formulas are the central result of this work. We note that in the course of our 
investigation, we discovered that the analog of the result f)2.22p has been found in [IT] 
in the setting of abstract graph theory and using the formalism of Grassmann variables. 
However, the proof is based on the notion of generalized star product [8] and requires 
a rather involved proof by induction on the size of the graph. Here, it is obtained 
directly by simple linear algebra and ready to use for computations (either analytical or 
numerical) . 
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2.2 Discussion 



We have checked that our formula reproduces known results obtained by other methods 
for simple graphs (star-triangle, etc.) [TJ [8j [10]. In the following, we present in detail 
the computation for new graphs, especially in 3D, for which the previous methods are 
impractical. Our method presents several advantages compared to previous ones. First, 
as just mentioned, it is computationally easier and one does not have to worry about the 
sequence of steps used in iterative methods where one has to make sure that fusing two 
given vertices and then a third gives the same results a fusing the first and third and 
then the second (cf [10]). The only task involved is the inversion of a matrix and there 
are well-known efficient methods both analytically and numerically. Then, we have an 
explicit formula which shows the location of the poles of S to t (on top of the possible ones 
from the local matrices which are given data in our approach). They are solutions of 

det(E(p) -S {22 \p)) =0. (2.24) 

This is important as these poles play a fundamental role in the computation of physi- 
cal quantities like the conductance in quantum systems defined on graphs (see [5j [TO]). 
Finally, for quantum systems on compact graphs, i.e. without external edges, the same 
equation provides the allowed modes on the graph. In this respect, f !2.24p is the gener- 
alization to an arbitrary compact quantum graph of the of the quantization equation 

e 2ikL = 1 , (2.25) 

for a particle in a box of length L. The matrix S^ 22 \p) accounts here for the one particle 
scattering occurring at the vertices. In the theory of integrable systems, this type of 
equations is sometimes called Bethe ansatz equations. However, here we emphasize that 
it is not related to such an ansatz. In condensed matter physics, the information given 
by this equation together with the dispersion relation of the model provides the basis of 
band structure analysis. 



2.3 Properties 

To be consistent, our general formula should not depend on the numbering of the inter- 
nal edges or vertices (internal permutation) and should transform appropriately under 
a permutation of the external modes (external permutation). Let II be an external per- 
mutation acting on A(p) and P an internal permutation acting on B(p). It is easy to 
see that this induces the transformations 

sM(p) ->ns , ( 11) (p)n- 1 (2.26) 

S {12) (p) ^US {12) ( P )P- 1 (2.27) 

S {21) (p) ^ PS^ip)^ 1 (2.28) 

S (22 \p) ^PS {22 \p)P- 1 (2.29) 

E(p) -> PE(p)P~ 1 , (2.30) 

producing Stat —> nS , to tiT~ 1 as it should. Therefore, in examples or applications, one 
can always fix a convenient numbering of edges and vertices and work up to an external 
permutation. 



6 



In view of physical application, we must also be concerned with the properties of 
S to t- We have seen already that S a {p)S a (—p) = l Na . This implies 

Stot(p)S tot (-p) = l Ne . (2.31) 

To see this, note that the block matrix made of (12.141) - (12. 171) is related to S(p) given in 
§2JM by 

S(p) 



{s^(p)\ s^(p) )- ps[p)P (2 - 32) 



where P is the permutation matrix defined by 

PA(P) = ( ) . (2 .33) 

Then by direct calculation and upon using S(p)S(—p) = 1^+2^ and E(p)E(—p) = JL 2 Ni 
we get 

StMSU-p) = l Ne +S (12 \p) (E(p) - S^ip))- 1 M (E(-p) - S^i-p))' 1 S^ 21 \-p) 

(2.34) 

where 

M = l 2Ni - (E(p) - sW(p)) (E(-p) - S^(-p)) - (E(p) - sM(p)) S^(-p) 
-S<&(p) (E(-p) - S^i-p)) - sW(p)SW(-p) 
= 0. (2.35) 

Now the local scattering matrices can be required to have additional properties, like 
unitarity. This is the case in particular if they arise from non-dissipative local boundary 
conditions emerging from self-adjoint extensions of the free one-dimensional Hamiltonian 
(see e.g. [7]). One then has unitarity 

SUp) = S-\p). (2.36) 

Following the same type of argument as above, one finds that Stot is a l so unitary. 

We finish this section by providing a few properties of E(p). It is symmetric and we 
have already seen that E(p)E(—p) = fl^- 111 particular E(0) 2 = ~S. 2 Ni so its eigenvalues 
are ±1 and are equally degenerate. Also, E(0) is a permutation matrix and E(p) is 

a generalized permutation matrix (with coefficients of the type e ip * ) which can be 
written as a product of a permutation matrix and a diagonal matrix 



where 



with 



E(p) = D(p)E(0) = E(0)D(p) , (2.37) 

N N a/3 
a,/3=l j=l 



D{p)D{q) = D(p + q) , p, q E C . (2.39) 
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3 Including loops 

The case of loops attached to single vertices can be treated with minor modifications in 
our formalism. Essentially, the idea is again to see a loop attached to a given vertex a 
as arising from the gluing of two edges attached to this vertex. This will be most easily 
incorporated in the general formalism if we use the following trick for notations. Let 
N a a 7^ be the number of loops attached to vertex a. To each loop j, j = 1, . . . , N aa 
correspond two modest a^f-iip) an d a< 2j{p) which are related by 

a~-i(p) = e- ipd r a ™(-p) , j = l,...,N aa . (3.1) 
We collect these modes in two-component vectors 

of (P) = ( %\p) ) > j=l,-,Na*. (3.2) 

We denote all the components of the local scattering matrix S a (p) related to the loop 
modes by s°f jk (p), j = 1, . . . , 2N aa} k = 1, . . . , N af3 ; s^ jk (p), j = 1, . . . , N al3 , k = 
l,...,2N aa ] and s%?j k (p), j,k = 1, . . . , 2N aa . Mimicking (I3.2p . we then define, for 

«S*(p) = f ^t 1 '^ J > j = l,--;N aa ,k = l,...,N a0 , (3.3) 

\ S a;2j,k\P) J 

= ( s Z,2^i(p) «S,a(p) ) , j = 1, . . . , A^ a/ 3, k = 1, • • • , N aa , (3.4) 



and also, 



3 a;jk 

Finally, we define 



( n \ _ ( S a%j-l,2k-l(P) S a°2j-l,2k(p) \ ■ u _ 1 /V fq ^ 

\ S a:2j,2k-l\P) S a;2j,2k\P) J 



e ipd f , if (3 ^ a and N af3 ^ 0. 



e f (p) = < . / o 1 \ (3.6) 
V^" [ i o ) , if/5 = aan d TWO. 

With all this, the relations defining scattering and propagation on the graph take the 
same form as before (cf (12. ip and (12.21) ) 

N N al 

<(P) = E E «S*(P) ^("P) Vj = 1, . . . , iV a , ; V/5 = 0, 1, . . . , iV (3.7) 

7 =0 fc=l 

and 

af (p) = ef(-p) af (-p) Vj = 1, . . . , N af3 ; V/3 = 0, 1, . . . , N . (3.8) 



5 Again, the choice of numbering is for convenience only and is irrelevant to the final results. 



Therefore, all the formalism and the results developed in section 12.11 hold in the same 
form, provided one substitutes e" /3 (— p) for e~ ipd j in (I2.12p . One should not be de- 
ceived by the apparent similarity of the results with or without loops. In general, the 
consequences of adding a loop in a given graph can be drastic. 

However, as the formalism suggests, allowing for loops in graphs opens the possibility 
that two topologically completely different graphs can have exactly the same total scat- 
tering matrix. This is illustrated on the example below. In particular, this shows that 
the uniqueness of the inverse scattering problem, as discussed in [9] , does not extend to 
the case of graphs with loops^|. 

We consider the two graphs depicted in Figure [T] below. To illustrate our notations, 
we have displayed the modes involved in the construction, dropping the p-dependence 
for conciseness. They are topologically completely different, one being a triangle with 
one external edge attached to each vertex and the other being a single vertex star graph 
with three external edges and three loops attached to it. Note that for the triangle, we 
drop the unnecessary Latin subcripts since N a/3 = 1 for all a = 1, 2, 3 and (3 — 0, 1, 2, 3, 
(3^ a. 




Figure 1: Two topologically different graphs with the same total scattering matrix. 
Left: triangle. Right: Star graph with loops. 



We assume that the scattering and propagation data is given as follows (we drop 
again the p dependence for clarity), 



For the triangle, 




4 Uniqucness is only guaranteed if one requires in addition that the number of vertices is maximal 
(cf Theorem 4.6 in 0). 



9 



giving the four blocks as defined in ( I2.14l) - (12.17l) in the form 



5 (H) 



c 00 

s l 











S 2 











c 00 
S 3 



5-(22) 



/ 


-22 
S l 


-23 

s l 














\ 




-32 
S l 


c 33 
























a 11 
S 2 


S 2 


















-31 
S 2 


c 33 
S 2 
























a 11 


S 3 




V 














c 21 
S 3 


-22 
5 3 


/ 



S (12) 



c 02 


c 03 




















-01 
S 2 


c 03 
S 2 




















-01 
5 3 


-02 
S 3 



S (21) 



( 


-20 
S l 








\ 




c 30 















-10 
S 2 












c 30 
S 2 















-10 
S 3 




\ 








-20 
S 3 


/ 



(3.11) 



and 



E, 



/ 








e -ipd 12 











\ 
















e -ipd 13 









e -ipd 12 




































g— ipd 23 









g— ipd 13 
















\ 











g— ipd 23 








/ 



(3.12) 



For the star graph, 



/ 



T 



t oo 

l l;ll 




xlO 
l l;ll 

t 10 

l l;21 

l l:31 



i00 
l l:22 



t 10 
l l;12 
-10 
L l;22 

t io 

l l:32 






xOO 
l l;33 

t io 

l l;13 
£10 

h;23 

t io 

l l;33 



.01 

l l;ll 

.01 

*1;21 

.01 

H;31 

t 11 

H;ll 

h;21 
*1;31 



.01 

U l;12 

.01 

l l;22 

.01 

l l;32 

t 11 
l l;12 

t 11 
U l;22 

l l;32 



.01 

h:13 

.01 

H;23 

.01 

H;33 

*1;13 

h;23 

t 11 

H;33 



T (ll) 


T (12) 


T (21) 


T (22) 



(3.13) 



and 



( 



\ 




















e ~ipd\ l 


































































g— ipd 
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(3.14) 



The lengths of the internal edges of the triangle are related to the lengths of the loop 
in the star graph by 



d 



12 



di 



ii 



d 



23 



4 1 



d 



23 



4\ 



(3.15) 
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00 
1;11 

01 
1;11 

01 
1;21 

01 
1;31 

,10 
4:11 



,10 
4:21 



,10 
4:31 



4:11 



t 11 

4:21 



c 00 ,00 
' 4;22 



_00 ,00 
^2 ) 4:33 



.00 



and the following relations for the scattering data hold, showing in particular the matrix 
structure defined in (13.31) - (13.51) in the case of loops, 

(3.16) 
(3.17) 
(3.18) 
(3.19) 



(«? 0) ) fg a =(0 o),t° 1 ] 13 =( 5 f 0) 
(0 ) , t% 2 ={ 5 f ) , t? 23 = ( o ) 

(00), t?* 3 2 = ( Bf ) , t' 



c 20 

s l 



c 30 

s i 



.22 



1;32 

10 
) 1:12 



01 
1:33 



t 



.10 



,10 

4; 13 



,10 

4;22 



,10 

4:32 



.30 

2 

] 






,10 

4; 23 



,10 

4:33 



.11 



t 11 

4; 12 



.13 



q31 
S 2 



t 11 

4:22 



c 33 
5 2 








.22 



(o 






c 20 
5 3 



c 10 
S 3 

t 11 

4; 13 



) 4;23 



.01 



.23 



t 



11 

1:31 



.32 







t 



11 

1;32 



.12 



t 



11 
1:33 







_33 
5 1 










.21 
'3 



.11 



(3.20) 
(3.21) 
(3.22) 
(3.23) 
(3.24) 
(3.25) 



The fact that these two graphs give rise to the same total scattering matrix follows from 
the fact that their scattering data are related by an internal permutation 



(3.26) 



/ 1 














o\ 








1 




















1 























1 





1 














\o 











1 


0/ 



such that 

p g(22) p-1 

Then, 



T (22) 



PE+P- 



E x 



£(12) p-1 _ T (12) ^ p £,(21) _ p(2l) _ (327) 



ntriangle 
tot 



,S(ii) + S (i2) p-^PEtP-'-P P- 1 )' 1 P 
T (u) + T (i2) _ T (22)y 1 T (2i) 



S, 



star 
tot 



(3.28) 



4 Platonic solids 

In this section, we illustrate the freedom on numbering and the use of formula ( 12.221) on 
the convex regular polyhedra known as Platonic solids (tetrahedron, cube, octahedron, 
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dodecahedron, icosahedron) [16] . Once the scattering matrix is known, physical quanti- 
ties associated to the graph can easily be computed, such as the conductance, using the 
formalism developped in [3] . The calculation essentially relies on the pole structure and 
the general techniques have been explicited in [TO] . 

We carry out explicit calculations in the case of the tetrahedron and the cube. This 
choice is primarily motivated by aesthetic and academic criteria rather than any par- 
ticular practical application. It also shows the computational advantage of our method 
over recursive ones on rather involved graphs. More precisely, we consider graphs whose 
internal edges and vertices correspond to Platonic solids and for which exactly one ex- 
ternal edge is attached to each vertex. This corresponds to N a p = 1, a = 1, . . . ,N, 
(3 = 0, . . . , N, a (3. Note that the condition of regularity yields d± = d for all 
a, j3 = 1, . . . , N. Also, all the vertices are connected to the same number of vertices so 
v a = v is the same for all a — 1, . . . , N. N is even for all those graphs. Finally, from the 
general theory of graph colouring, see e.g. [17J, it is known that we can assign a label (or 
colour) a G {1, ...,;/} to the edges connected to the same vertex in a way compatible 
with the graph i.e. , in colour terms, such that no two edges connected to the same 
vertex have the same colour and each edge can only have one colour. This allows us to 
define functions n a , a G {1, . . . , N} from {1, . . . , z/} to {1, . . . , A^} such that n a (a) = (3 
if and only if (3 is connected to a by the edge labelled a. We use the convention a = 
for external edges and set n Q (0) = for all a G {1, ... , A^}. By construction, we have 
the following properties 

n a (a) = (3 -v^ np(a) = a , n a (a) = np(a) <^ a = (3 , n a (a) = n a {b) a = b . (4.1) 

In view of formula (12.221) . the main object of interest is E(p) — S t - 22 \p) which we seek 
to invert. With our notations, we get 

N v N v 

a=l a=l a=l a,b=l 

where the local matrices read 

S a (p) = E a+ i, b+1 ® < 6 (p) , a = 1, . . . , N . (4.3) 

a,b=0 

For later convenience, we define a reduced scattering matrix containing only the infor- 
mation about scattering on the internal edges 

V 

S r :\v) = £ E a , b ® 5 ^(p) , a = 1, . . . , N . (4.4) 

a,b=l 

Let us also define 

N 

£ a = J2 E «,n a (a). (4.5) 
a=l 
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Then E{p) = e %v 2_^^ a ® ^ aa an< ^ f rom the general properties of E (or by direct 

a=l 

calculation) we find 

£ a = £ t a = £~ 1 , a = l,. ..,!/. (4.6) 

Therefore £ a is diagonalizable with eigenvalues ±1 each degenerate y times and with 
eigenvectors v £ a = 7^(e a + ee na („)), e = ±1, a < n a (a), forming an orthonormal basis. 

4.1 Tetrahedron 

For the tetrahedron, N = 4, v = 3 and the matrices £ a enjoy the additional property 

£ a £ b = £ b £ a , V a, 6= 1,2,3, (4.7) 

due to the fact that 




Figure 2: Tetrahedron with an example of numbering. 



From (14. 7p . they can be diagonalized simultaneously. As already explained, to fix ideas 
we can fix a numbering without loss of generality since we work up to permutations. In 
the present case, changing the edges and or vertices numbering amounts to interchanging 
the £aS. From the figure we obtain 

/ 1 \ / 1 \ 

0010 1000 
1_ 0100 ' 2 ~ 0001 

\ 1 / \0 1 0/ 



£* 



( 1 \ 

1 

10 

V 1 / 



(4.9) 
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and a diagonalizing matrix is 



1 

T =2 



-1 1 -1 1 \ 

1 -1-11 
-1-1 1 1 

Viiiiy 



1-1 



(4.10) 



So far, we haven't taken advantage of the geometry and its symmetries. The scattering 
can still be different from vertex to vertex (as labelled by the index a on the local 
matrices) and at a given vertex, the scattering from edge a to edge b needs not be the 
same as the scattering from edge a to edge c say. Clearly, this does not respect the natural 
symmetry of the underlying graph. One can impose that the local scattering matrices 
be the same for all vertices i.e. S a (p) = S(p) and in particular S Ted {p) = S red (p) for all 
a = 1, 2, 3, 4. This already greatly simplifies the problem of inversion. Let r = T ® 1 3 
and D a = T£ a T~ l then 

3 

r(E(p) - S^ip))^ 1 = e~ ipd J2 D a® E aa - 1I 4 ® S red (p) (4.11) 

a=l 

The matrix on the right-hand side is a block diagonal matrix made of four 3x3 blocks 
essentially determined by S red 

4 

r(E(p) - S^ 22 \p))t- 1 = ^E aa ® (e-^/ Q - S red (p)) , (4.12) 



with 






Is -(4.13) 



Thus, the problem is reduced to inverting 3x3 matrices. In particular, the poles of S to t 
are solutions of 

det (e~ ipd I a - S red (p)) =0 , a = 1, 2, 3, 4 . (4.14) 

We now turn to the explicit calculation of S to t m the case where the vertices are de- 
scribed by scale invariant local matrices (independent of p) capturing universal features 
of scattering. In our case, each local matrix is the same 4x4 scale invariant matrix 
whose explicit form has been classified in [3J. Note also that we can take further ad- 
vantage of the symmetries of the underlying geometry here by imposing for instance 
that the scattering be invariant under a rotation of 7r/3 around the axis passing through 
a vertex and the centre of the opposite face. Physically, this means that an incoming 
particle from the external edge of a vertex has the same probability of being transmitted 
to any one of the internal edges attached to this vertex. Mathematically, this amounts 
to requiring that S satisfies 

1 ° Vi M=*. MS) 



J ) \ J 
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where 



J 



1 

1 

1 



Putting everything together, we find two possible local scattering matrices 



/ _I I 

' 2 2 2 2 
1 111 



Si 



1 1 \ / _1 I 1 1 \ 

9 / 9 9 9 9 



2 2 2 

1 1 1 



2 2 

1 1 

2 2 



2 2 
1 1 



, S 2 



/ 



V 2 2 

In the first case, we compute the S to t{p) as 
1 



2 2 2 2 
15 11 



6 1 5 6 



I I _I _l 

\ 2 6 6 



5 

6 6 



sUp) 



(-2 (e- 3ipd + e- ipd - 1) I 4 + e -' ipd (e- ipd + 1)A) 



where dip) = (2e- 2ipd + e~ ipd + l)(2e- ipd - 1) and 



A 



( 1 1 1 \ 
10 11 
110 1 

Vi i i o ) 



The poles of this matrix are given by 

- ipd = x with xe\^, 



1 -l+iv 7 ^ -l-i-v/7 



In the second case, we obtain 
1 



(4.16) 



G 2 (p) 



where G 2 (p) = (6e~ 2ipd - e~ ipd - 3)(2e~ ipd - 1), leading to the poles 

1 1 + V7S 1-V73 



e ipd = x with x e |-, 



12 



12 



}■ 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(-2 (-6e~ 3ipd + Ae~ 2ipd + 10e"^ - 6) I 4 + 3e~ ipd (e~ ipd - 1)A) , 

(4.21) 



(4.22) 



4.2 cube 



For the cube, N = 8 and z/ = 3 and the matrices S a also commute. So one can perform 
the same analysis as before. 
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Figure 3: Cube with an example of numbering. 



Based on figure El we get explicitly 
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1 



















1 
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, £2 = 
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1 





























1 














1 































1 


























1 



















1 

























1 

















1 








0/ 



(4.23) 



S, 



/o 











1 








o\ 

















1 


























1 


























1 


1 


























1 


























1 

















\o 








1 











/ 



(4.24) 
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and a diagonalizing matrix is 



V 



2V2 



/ 
1 


1 
1 


1 

— 1 


1 
1 


1 

— 1 


1 

— 1 


1 
1 


1 

— 1 


1 \ 




1 


1 


-1 


-1 


-1 


-1 


1 


1 




-1 


1 


1 


-1 


1 


-1 


-1 


1 




-1 


-1 


-1 


-1 


1 


1 


1 


1 




-1 


1 


-1 


1 


-1 


1 


-1 


1 




-1 


-1 


1 


1 


-1 


-1 


1 


1 




1 


-1 


-1 


1 


1 


-1 


-1 


1 


V 


1 


1 


1 


1 


1 


1 


1 


1 / 



(V 



-l\t 



(4.25) 



Again assuming that the local scattering matrices are the same at all vertices, we get 



V(E(p) - (p^V- 1 = e^ d J2 A a® E aa - Is ® S red ( P ) 



(4.26) 



a=l 



where A a = V £ a V 1 and V = V ® H3. This is a block diagonal matrix and the problem 
is reduced to inverting 3x3 matrices, 



(E(p)-S^\p))- l = V- 1 



Y^E aa ® (e^ d X a - S red {p)Y l 



a=l 



V. 



(4.27) 



where 



Ji = i 3 , x 2 = 

1 
1 
0-1 



-1 









( 1 












1 










-1 













1 J 




V 





1 J 





-10 
0-10 
1 



-1 





^ 




f 1 










1 










-1 


• 








-1 J 




V 








,(4.28) 
H 3 • (4.29) 



We turn to the explicit computation of the total scattering matrix in the two cases (14.171) 
describing scale and rotation invariant local scattering at the vertices. In both cases, we 
find the following structures for Stot'- it is a linear combination of matrices in the abelian 
group generated by the £'s with coefficients being polynomials in e~ tpd . For j = 1,2, 

Stotip) = a o(p)^s + a{(p)£i + 4(p)^2 + 4(^)^3 + a{{p)£i£2 + a{{p)£i£ z + a J 6 (p)£ 2 £ 3 
+4{p)£ 1 £ 2 £ 3 . (4.30) 
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In the first case, we find 

g + e ~i P d _ 8e ~2i P d _ 5e ~3i P d _ 4Q e -4i P d + 4e -5i P d _ 32e -6ipd 

a iW - 4 (_! + e -2i P d + 8e -upd + i 6e -6 ip d) > t 4 - 32 ) 

= 4 _ 16e - 2 ^ . ( 4 -33) 

= 4 _ 16e ^ P d . ( 4 -34) 

a 4 (p) = - j^— ^ - - 2e _ 2ipfl! _ 4e „ 3 i p(i _|_ g e -4j P d) ' (4.35) 

1 e~ ipd (1 - 9e~ ipd + 2e~ 2ipd ) 

a s(p) - ~ _)_ e -jpd + 2 e - 2i P d - 4 e - 3i P d + 8e~ 4i P d ) ' (4.36) 

e~ ipd (1 + 9e~ ipd + 2e- 2ipd ) 

a eW _ 4 (-1 - e _i P d + 2e~ 2i P d + 4e~ 3i P d + 8e- 4ipd ) ' (4.37) 

g-ipd + 19 g-3ipd + 4e -5ipd 

- ~ 4 (_1 + e -2rpd + gg- 4l pd + 16e -6ipd) • ( - 38 ) 

The poles of the scattering matrix can be then computed. They are given by 

e-** = z with xG {±l,±ii^,±iz^}. (4.39) 
In the second case, we find 

72 + g e -ipd _ 44 Qg-2ipd _ 45e -3ipd + 72 gg-4ipd + ggg-Sipd _ 2 gg e -6ipd 

W) - 4 (_ 9 + 7 3e -2^ _ 184e -4ipd + iu e -^ d ) ' ( ' ^ 

3e" ipd (15 - 67e- 2ipd + 60e~ 4 ^) 

a iW = ~4 (-9 + 73e- 2i P d - 184e- 4 ^ + 144 e - 6 ^) ' ^ ^ 
3e~ ipd 

= 4 _ lQe -2,pd » (4-42) 

a 5<P) = 4 _ 16e - 2 ,p, , (4-43) 

3e~ ipd (-1 + 3e~ ipd + 2e~ 2ipd ) 

° 4 ^' = ~4 (3 - e-^ - l8e~ 2i P d + 4e" 3ipd + 24e~ Ai P d ) ' ^ ' ^ 

3e- ipd (-l + 3e-^ + 2e- 2ipd ) 



a 2 7 (p) 



4 (3 _ e -ipd _ ig e -2«pd + 4 e -3ipd _|_ 24e _4ipd ) ' 

f4 46) 

4 (3 + e -«pd _ ig e -2^ - 4 e -3ipd _|_ 24e _4ipd ) ' ' 

i i . (4.47) 

4 (9 - 73e- 2ipd + 184e~ 4ipd - 144e~ 6ipd ) 
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The poles of the scattering matrix are given by 

e-* = x with s G {±i,±!+^,±i^_^}. (44g) 

5 Conclusion and outlooks 

In this paper, we have presented a direct method for the computation of the total 
scattering matrix of an arbitrary finite noncompact connected graph given its topology, 
metric structure and local scattering data at each vertex. The method uses the formalism 
of quantum modes as our initial motivation was the study of quantum fields on graphs. 
This resulted in a simple and direct algebraic derivation of formula (I2.22p . We have 
also shown that the case of loops is easily incorporated in our method. This has been 
illustrated with an explicit example whose purpose was also to point out that the inverse 
scattering problem on graphs does not have a unique solution in general for graphs with 
loops. 

We want to stress that in the present paper, the modes as we called them, appear 
more as convenient labels than true quantum field theoretic objects. This has two 
consequences. First, our results are ready to use for applications in quantum field theory 
on graphs by simply promoting the modes to generators of the RT-algebra [14] . Second, 
it means that our results hold in complete generality for abstract graphs with or without 
loops. In this respect, the present results provide an extension of the results in [TT] to 
the case of loopcl 

Finally, this paper lays the ground to applications to transport problems on arbitrary 
graphs in the spirit of the studies performed in e.g. [31 El EH E] . Indeed, it provides one 
with the central ingredient which is the total scattering matrix together with its pole 
structure. Once this structure is known, the calculation of physical data such as the 
conductance between external edges is rather direct, see e.g. [10J: we will return to 
these questions in the near future. 
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